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ABSTRACT 



Context. The central field of the Andromeda galaxy (M 31) has been monitored, using the Chandra HRC-l deteclor (about 0.1- 10 keV 

energy range) from 2006 to 2012 with the main aim to detect X-rays from optical novae. We present a systematic analysis of all 

X-ray sources found in the 41 nova monitoring observations, along with 23 M 31 central field HRC-I observations available from the 

Chandra data archive starting in December 1999. 

Aims. Based on these observations, we studied the X-ray long-term variability of the source population and especially of X-ray 

binaries in M 3 1 . 

Methods. We created a catalogue of sources, detected in the 64 available observations, which add up to a total exposure of about I 

Ms. To study the variability, we developed a processing pipeline to derive long-term Chandra HRC-I light curves for each source over 

the 13 years of observations. In the merged images we also searched for extended X-ray sources. 

Results. We present a point-source catalogue, containing 318 X-ray sources with detailed long-term variability information. 28 of 

which are published for the first time. The spatial and temporal resolution of the catalogue allows us to classify 115 X-ray binary 

candidates showing high X-ray variability or even outbursts in addition to 14 globular cluster X-ray binary candidates showing no 

significant variability. The analysis may suggest, that outburst sources are less frequent in globular clusters than in the field of M 31. 

We detected 7 supernova remnants, one of which is a new candidate and in addition resolved the first X-rays from a known radio 

supernova remnant. Besides 33 known optical nova/X-ray source correlations, we also discovered one previously unknown super-soft 

X-ray outburst and several new nova candidates. 

Conclusions. The catalogue contains a large sample of detailed long-term X-ray light curves in the M 3 1 central field, which helps to 

understand the X-ray population of our neighbouring spiral galaxy M 3 1 . 

Key words. Galaxies: individual: M 3 1 - X-rays: galaxies - X-rays: binaries - Galaxies: stellar content - Catalogs 



1. Introduction 

M 3 1 is an ideal target for X-ray observations since it is the clos- 
est spiral galaxy and similar to our own galaxy. The outside view 
and the distance of M31 enable us to monitor several hundred 
X-ray sources in just one telescope pointing. The Einstein X-ray 
observatory ( e.g. Trinchieri & Fabbia no_1991) an d the ROSAT 
satellite (Pri mini et al.||1993 ; Supper etal ||1997| |200r) gave a 
first insight into the X-ray population of M31, followed by the 
instruments of Chandra (e.g. Kong et al.ll2002b ; K aaret|[2002 
|Wimamsetal .'2004'; Voss & Gilfanov 2007a) and XMM-Newton 
(e.g. |Pietsch~et al. 2005b; Stiele et al.^2008, 2011,) with better 



been studied. The main class of M 3 1 X-ray sources are X-ray 
binaries (XRB), in which a compact object like a white dwarf, 
a neutron star or a black hole accretes matter from a compan- 
ion star Such sources can show highly variable or transient be- 
haviour (e.g. Osborne et al. 2001 , Williams et al. 2006 ). The sub- 
class of globular cluster X-ray sources (GIC), which are mostly 
low-mass XRBs (LMXBs) is covered in several papers (e.g. Di 
Stefano et al.|2002[[lYudolyubov & Priedhorsk y 2004; P eacock 



spatial resolution and deeper observations. 
Several individual X-ray sources in M31 have been discussed 
in previous publications (e.g. sources showing X-ray dips, pul- 
sations, bursts or ultra-luminous X-ray outbursts, for references 
see sectionHland Table 1). In addition many source classes have 

* The primary source catalogue table (Table 1) and the variability 
table (Table 2) will only be available in electronic form at the CDS 
(|http://cds.u-strasbg.fr/^. 



etal . 2009; Barnard et al. 2012b). A further subclass of XRBs 
are defined by their super-soft X-ray spectrum (SSS) (e.g. ffl 
IStefano et al.|2004| i. Most of t hese SSSs have been found to be 
associated with optical novae ( Pietsch et al.||2005a 2007 Orio 
I 2006| ). Another class of X-ray sources are supernova remnants 
(SNR), which can be spatially resolved in M 3 1 using X-ray, op- 
tical or radio observations ( ,Kong et al., 2002a Sasaki et al. 2012 
and references therein). 

Beginning in 1999 the local group galaxy M31 (Andromeda 
Galaxy) has been observed regularly with the Chandra X-ray 
Observatory. The satellite High Resolution Mirror Assembly 
(HRMA, see Weisskopf 1988|l, was frequently pointed at the 
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Fig. 1. Optical image of the Andromeda Galaxy (DSS2 blue) 
with X-ray source regions from the Chandra HRC-I catalogue 
(~ 50% encircled count fraction radius) and a circle marking the 
outer rim of the detector FOV (image obtained from the NASA 
Sky View website). 



centre of M31, which allows us to examine the long term 
evolution of the intensity of sources within the central field 
of the galaxy. We used 64 observations of the Chandra High 
Resolution Camera in imaging mode (HRC-I, see Murray et al. 
[T997| ) from 1999 to 2012. The single Chandra HRC-I obser- 
vation with the longest exposure time is discussed by Kaaret 
( |2002] l 



Williams et al. (200411 discuss snapshots, covering the 



whole galaxy until 2002. 
The observations from 2006 to 2012 are part of an extensive 
monitoring program (41 observations with about 20 ks expo- 
sure), to investigate X-ray sources associated with optical novae. 
These observations have been searched for novae by Hen zeetaT] 
( 2010| l and Henze et al7| ( ,201 1) . We now use these monitoring ob- 
servations, along with all previous archival HRC-I observations 
to systematically analyse the time variability of all sources, de- 
tected in the M 3 1 central field with Chandra HRC-I. 
We used the HRC-I detector, which has a wide field of view (32' 
X 32') and a very high angular resolution (below 1") for sources 
close to the pointing of the telescopaH The dense HRC-I moni- 
toring of M 3 1 allows us to create the deepest (about 1 Ms of total 
exposure) catalogue of X-ray sources at the centre of M 3 1 and 
provide a detailed long-term variability study for each source. 
Figure [T] shows an image of the digitized sky survey 2 (DSS2) 
taken with a blue filter and the approximate HRC-I field of view 
(FOV). 

2. Observations and data reduction 

As mentioned earlier, there are 64 Chandra HRC-I observa- 
tions to date, pointed to within 3' from the cent re of M31 (RA: 
00^ 42™ 44.33\ Dec: -h41° 16' 07.50", J2000, ISkrutskie et al. 



|2006[ ). We excluded two ~1 ks observations with high back- 
ground count rates (ObsID 00270 and 01569), which would not 
have contributed new sources to the catalogue and would have 
deteriorated the signal to noise ratio, when merging the individ- 
ual observations. Table IaTI contains information on all HRC-I 
observations used. The observations from 1999 to 2002 mostly 
have an exposure of ~1.2 ks (one with 2.5 ks and one with 
5.2 ks), except for ObsID 1912 (46.7 ks), which was used by 
Kaaret (2002) for his investigation. From 2004-12-06 till 2005- 
02-21 there were six HRC-I observations of the M31 centre, 
with exposures between ~20 and ~50 ks to monitor the vari- 
ability of M31*, the source at the M31 nucleus. The period 
from 2006 to 2012 is covered by the ~20 ks HRC-I observa- 
tions of the XMM-NewtonI Chandra nova monitoring program 
and additional Chandra guaranteed time observations of the 
Max-Planck-Institut fiir extraterrestrische Physik (MPE). In this 
period we have rather regularly spaced data points with about 
the same exposure. This monitoring delivers the best input for 
the time variability studies. 

2. 1 . Source detection 

The data for each observation were reprocessed using the 
Chandra Interactive Analysis of Observations software pack- 
age (CIAO, [Fruscione et al.|2006|l version 4.4 an d the Chandra 
Calibration Database (CalDB, |Graessle et al.| |2007[ l version 
4.4.7, to receive a level 2 event list according to the Chandra 
Standard Data Processing (SDP). For each observation we cre- 
ated a light curve of the whole detector and removed events at 
times of high background from the event file (good time interval 
(GTI) filtering) using the Munich Image Data Analysis System 
(MIDAS) software package provided by the European Southern 
Observatory (ESO). We only used events in the pulse invariant 
(PI) channel range from 48-293, which according to the CIAO 
science threads increases the contrast between sources and back- 
grouncj^ 

As the point spread function (PSF) of the Chandra HRMA 
causes higher source-blurring with increasing off-axis angle, we 
used five images with different resolution (similar to Kaaret] 
2002 1. This reduces computation time, since the detection does 
not have to run on the entire detector field of view at high- 
est resolution. The highest resolution image covers 200"x200" 
with 0.13" pixel size (HRC-I detector resolution). The sec- 
ond, third and fourth image cover 400"x400", 800"x800" and 
1400"xl400", with abinning of 2x2, 4x4 and 8x8 detector pix- 
els. The fifth image covers the entire detector and has a 16x16 
binning. 

The size of the PSF increases by a factor of ~ 30 from the centre 
to an off-axis angle of ~ 15'. where the effective exposure de- 
creases to about 70%. This means that the spatial resolution and 
the sensitivity of the catalogue depends on the off-axis angle of 
a source. 

Using CIAO tools we created an exposure map for each image, 
correcting for vignetting effects of the HRMA and for the ef- 
ficiency of the detector. After testing various background mod- 
els (by subtracting background maps from images and inspect- 
ing residuals) we found that the XMM-Newton Science Analysis 
System (SAS, [Gabriel et al.||2004[ ) task espliiiemap delivered 
the best results for modelling background contributions and the 
diffuse emission of M31. To separate point-source contribu- 
tions, we did a detection run with the CIAO tool wavdetect 
with the sensitivity parameter sigthresh - 10""*. We chose a 



http://cxc.harvard.edu/proposer/POG/html/HRC.html 



http://cxc.harvard.edu/ciao/threads/hrci_bg_spectra/ 
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slightly higher sensitivity than recommended in the CIAO sci- 
ence threads^ to avoid missing any detections. Later we removed 
detections below a significance of 3cr. To avoid spurious detec- 
tions at the detector edges, we used a detection mask, to exclude 
areas where the effective exposure is below 50% of the maxi- 
mum exposure, due to vignetting effects of the mirrors and de- 
tector efficiency. The scales of the wavelets were optimized to 
detect point sources in the relevant region of each image. The 
esplinemap tool used the wavdetect output to remove sources 
above a likelihood of 5cr from the image and then created a back- 
ground map, fitted to the re-binned image with 16 spline nodes. 
With these background maps as input we then ran wavdetect 
again on every image, with the same parameters as in the first 
run, to create a final source list for each image. 
We merged the output source lists of the five images. The lists 
cover regions from radius 0-95", 90-200", 190-400", 380-700" 
and 665-1315". Sources that were found in two images in an 
overlap-region of two lists were taken from the higher resolution 
image. 

Based on the wavdetect count rates, we calculated the 0.2- 
10 keV flux, correct ed for Galactic foreground absorption 
(Nh = 6.6 ■ 10^° cm^^, |Starketal.|l992| l towai-d M 31 and used a 
power-law spectrum with a photon index of F = 1.7. For the lu- 
minosity calculation, we ass umed a distance of 780 kpc to M 3 1 
( Stanek & Garnavich||1998| l and isotropic radiation, which has 
been used in many previous X-ray studies. Although there are 
more recent papers on the distance to M31 (e.g. Riess et al. 
20T2I 752 + 27 kpc). 



To improve the astrometry of individual observations we com- 
pared the source-list positions with the'Kaaret'(2002) catalogue, 
which is based on the longest individual observation and then 
corrected systematic offsets. The Kaaret (2002i catalogue has 
been astrometrically corrected after cross correlation with the 
Two Micron All-Sky Survey catalogue (2MASS, see Skrutskie 
|et al .'2006^). After this correction, the mean angular distance to 
the |Kaaret,(2002j sources was 0.23". Since Kaai'et,(2002) states, 
that systematic uncertainties are about 0.15", we estimated a 
conservative systematic Icr uncertainty of 0.38" for the cata- 
logue. 

2.2. Merged image 

We merged the individual observations to find possible faint 
sources, that have not been detected in the individual observa- 
tions. 

We corrected for the systematic offset against the Kaaret ( 2002| l 
catalogue by re-projecting the events for each observation. We 
again produced images at the five resolutions (see |2.1| i, but 
aligned to the coordinates of the centre of M 3 1 instead of the 
centre of the detector. Thereby, we ensure that a given pixel 
coordinate in one image coiTesponds to the same RA/Dec and 
pixel coordinate in the images of each observation. Therefore, 
we were able to add up the images pixel by pixel, instead of 
merging the event lists, which would have become too large to 
handle for some of the programs and would have used up too 
much computational resources. 

The images and exposure maps of the individual observations 
were merged. We created new background images, accepting 
only merged areas, for which the effective exposure is above 
70% of the maximum. This was necessary, because the corners 
of the detector fields do not overlap further off-axis in the 
merged image, as the roll angle of the satellite changed during 
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Fig. 2. Colour image of the merged HRC-I exposure maps, 
showing the entire detector FOV. The dashed black circle marks 
the outer limit of the catalogue (see Fig. [TJ. The dashed white 
circle includes the approximate area with > 70% effective expo- 
sure. The colour bar shows the effective exposure time value (in 
seconds). The innermost black region marks the approximate 1 
Ms effective exposure region. 



the course of the observations, leading to a fast drop in exposure 

(see Fig.|2]). 

For these images, we ran the detection pipeline, described in 

section IZTI 

The maximum exposure time of the merged image of all 62 

observations, we used, adds up to 1.02 Ms in the central area 



(see Fig. 14 15 16 1 



2.3. Spectral information 

In order to distinguish between XRBs with harder spectrum or 
SSSs with super-soft X-ray spectrum, we checked the hardness 
ratios (HR) of correlating X-ray sources from previous cata- 
logues (see section 3.2 1. In some cases we also checked hard- 
ness information from previous studies on individual sources, 
in which data from the XMM-Newton European photon imag- 
ing camera (EPIC, [Striider et aL][200T] [Turner et al.,,2003) de- 



tectors, the Swift X-ray telescope (XRT, Gehrels et al. 2004|i 



or the Chandra Advanced CCD Imaging Spectrometer (ACIS, 
Garmire et al.|2003| ) were analysed. Where we could not obtain 
spectral information from the identification with known X-ray 
sources, we extracted hardness information from one adjacent 
detection with Chandra ACIS, obtained from the Chandra data 
archivqj although there were many more observations available 
in most cases. We used observations with the ACIS-I or ACIS-S 
CCD array. 
We extracted ACIS source counts in the 0.2 - 1 keV and 1 - 



http://cxc.harvard.edu/ciao/tlireads/wavdetect/ 



http://cxc.harvard.edu/cda/ 
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Table 3. Photon count distribution observed with Chandra ACIS in the 0.2 keV to 10 keV photon energy range. 



No." RA(J2000) Dec(J2000) Observation' 



Photon counts (ct) 
0.2keV < ct < IkeV 



IkeV < ct < lOkeV 



47 


00:42:21.29 


+41 


15:52.6 


ACIS-I,ObsID 11275, 


2009-11-11 


65 


00:42:29.08 


+41 


15:47.6 


ACIS-I, ObsID 9523, 


2008-09-01 


127 


00:42:41.46 


+41 


16:17.9 


ACIS-I, ObsID 7064, 


2006-12-04 


128 


00:42:41.61 


+41 


14:37.2 


ACIS-I, ObsID 7064, 


2006-12-04 


154 


00:42:43.80 


+41 


16:12.7 


ACIS-I, ObsID 9524, 


2008-10-13 


159 


00:42:43.93 


+41 


16:10.8 


ACIS-I, ObsID 9524, 


2008-10-13 


161 


00:42:43.99 


+41 


16:37.1 


ACIS-I, ObsID 12974, 


2011-09-28 


162 


00:42:44.12 


+41 


16:04.3 


ACIS-I, ObsID 7138, 


2006-06-09 


180 


00:42:45.10 


+41 


15:42.8 


ACIS-S, ObsID 14197, 


2011-09-01 


189 


00:42:46.42 


+41 


16:10.2 


ACIS-I, ObsID 11279, 


2010-03-05 


222 


00:42:51.79 


+41 


17:27.3 


ACIS-I, ObsID 7064, 


2006-12-04 


244 


00:42:55.99 


+41 


17:21.5 


ACIS-I, ObsID 12161, 


2010-11-16 


247 


00:42:56.97 


+41 


20:05.1 


ACIS-I, ObsID 12160, 


2010-10-19 


259 


00:43:01.15 


+41 


13:17.7 


ACIS-I, ObsID 7064, 


2006-12-04 



57 
14 

271 
87 
28 
19 
23 
18 

108 
41 

141 

6 

48 

12 



50 
21 

413 

190 

42 

20 

15 

9 

151 
62 

120 

4 

49 

22 



Notes. *"* Source number in the catalogue, 
date). 



Observation, from which the spectrum was estimated (instrument, observation number (ObsID), 



10 keV band (see Table [3]l. The counts were extracted from re- 
gions with ~ 95% encircled count fraction PSF size and were not 
corrected for background contributions. For the ACIS-I observa- 
tions the background was negligible. In the ACIS-S observation 
(ObsID 14197, ~ 37.5 ks exposure) we extracted 9 counts in the 
background. The obtained count ratio given in Table [3] can be 
used to clearly distinguish hard X-ray spectra from super-soft 
spectra. In all cases the source counts in the bands were of sim- 
ilar order, indicating a hard spectrum. For a SSS spectrum most 
of the counts would be expected in the 0.2 - 1 keV band. SSS 
selection criteria have been discussed by Kong et al. ( 2002b| l for 
Chandra ACIS and by jPietsch et aL| ^^^ffoFXMM -Newton 
EPIC. 

For two outburst sources (No. 1 16 and No. 300, see section [4~T| , 
which were only observed with Chandra HRC-I, we give a rough 
estimate of hardness based on HRC-1 HRs, as has been success- 



fully demonstrated by Pietsch et al. (2007 ). 
The Chandra observatory proposers guida^ suggests this proce- 
dure. 



3. Analysis and results 

3.1. Catalogue compilation 

To produce the source catalogue from the detection lists of the 
individual and merged observations (for sources with signifi- 
cance above 3cr), we matched positions of sources in each ob- 
servation, using a script implemented in Pythorrl We consecu- 
tively merged the detection lists of the individual and merged 
observations, by adding sources, which were not correlating 
with any source position in the preliminary merged list within 
a radius of two times the 50% encircled count fraction PSF 
size. This conservatively exceeds the 3cr positional errors given 
by wavdetect and in addition accounts for systematic effects, 
which are caused mainly by the strong degradation of the PSF 
size with increasing off-axis angle. While building this list we 
excluded sources, which were only detected in one observation 
and had a source significance below 5cr. The parameters adopted 

^ http://cxc.harvard.edU/proposer/POG/html/chap7.html#tth_chAp7 
* http://www.python.org/ 



for each source are derived from the most significant detection 
with wavdetect. We then removed sources with a maximum 
source significance below Act and screened the list. During the 
screening, we removed one duplicate of a source, spurious de- 
tections around very bright sources, faint and far off-axis sources 
with bad positions and false sources found at the detector edges 
with unrealistic wavelet scales. Overall there were 11369 indi- 
vidual detections with wavdetect. After the merge and after 
applying the > 3cr cut, there were 384 sources in the catalogue. 
The further significance cuts and the screening reduced the num- 
ber of sources from 384 to the final 318. 

As the detection was optimised to find point sources, we also 
removed known supernova remnants (SNRs included in Sasaki 
|etal. 12012 1 from the catalogue, if the source was split into mul- 
tiple spurious detections, due to the extent of the SNR. We will 



discuss them separately in section 4.2 The wavdetect tool did 
not find the correct position for the far off axis neighbouring 
sources No. 20, 28 in the catalogue (with the merging param- 
eters), so we determined the detection properties from an ob- 
servation (HRC-1 ObsID 5925) where one of the sources was 
too faint to be detected. The properties of No. 18 and 21 were 
also determined from single detections (HRC-1 ObsID 5926 and 
5927) since they were too close together to be resolved during 
the automatic compilation of the catalogue. We also added the 
positions for the double nucleus of the galaxy (No. 164 and 165 
see Fig. [3]), from ObsID 10684, where the nucleus is resolved. 
The final catalogue includes 318 sources, 20 of which are only 
detected in the merged images. It is sorted ascending in right 
ascension (J2000). From this point onwards we shall refer to 
specific sources in the catalogue with their source number (see 
SRCJD column in section [33] ) and "No." in front. 
It should be noted, that the systematic positional errors increase, 
the further off-axis a source is located as the PSF of the Chandra 
HRMA degrades with increasing off-axis angle. Figures 14 15] 



16 show the source regions of the HRC-1 catalogue plotted 



over the merged images at different resolutions. The sources are 
marked by circles with ~ 90% encircled count fraction PSF ra- 
dius for 1 keV photons at the source position. 
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44.5 0:42:44.3 



44.1 



Fig. 3. Logarithmically scaled contour image resolving the M 3 1 
nucleus (ObsID 10684). Coordinate system in RA and Dec. 
Image smoothed with Gaussian function of 2 pixel FWHM and 
contours are chosen at 0.5, 1.1 and 1.5 counts per pixel to il- 
lustrate the separation of the two sources at the core of M 3 1 
(double core east No. 165 and west No. 164). No. 162, 169, 170, 
172 were not active in this observation. The region size is the 
90% encircled count fraction PSF. 



3.2. Correlations 

The catalogue was correlated with the X-ray catalogues, listed 
below (abbreviation, reference and assumed 68% confidence 
errors) to identify previous detections and classifications of 
sources: 



- KGP2002: Kong et al. (' 2002b| l, individual eiTors given 

- Ka2002: Kaaret, (2002) 0.35'' 



PFH2005: 
SHL2001: 



Pietsch et al. 



Supper et al. 



! 2005b I, individual errors given 
2001), individual errors given 



SHP97: [Supper et al.|(|19 97), individual errors given 
(2008), individual eiTors given 
(2011), indiv idual eiTors given 



SPH2008: Stiele et al 
SPH201 1 :, Stiele etal.. 
- VG2007: |Voss & Gilfanov| ( |2007a| l, 0.4" error 

K 0.5" error 
I, 0.5" error 



- WGK2004?Williams et al. 



- WNG2006: Williamset al. 

- HPH2011:He nzeetal 

- PHS2007 



E004 



. 2006 
( |2011|) 0.4" error 



jPietsch et al. ( 2007 1, individual errors given 

- Chandra Source Catalog (CSC): |Ev"ans et al.| ( [2010| ), individ- 
ual erro rs given 

- D2002: p^tefanoetal.|([2002l), 1" error 

- DKG2 004nDi Stefano et al.| ( |2004| ), 1" error 

- O2006: |OrJoR2006|l, 2" er ror 



- OBT2001: 

- SBK2009: 



Osborne et aLl(|200T]l I" er ror 

Shaw Greening et al. (2009|), 3.1" error 



- TF9 1 : jTrinchieri & Fabbiano ( 1 9 9l|rT.5" er ror 



- TP2004:'Trudolyubov & Priedho rsky|p004[ ), 1.5' 

- TP2006: |Trudolyubov et al.| ( |2006| 7l" error 



We considered two source positions as coiTelated, if their angu- 
lar distance (d) from each other meets the criterion 



d < 3.439 ■ Vo-i^ + o-2^ 



(1) 



where 0-^2 are the 68% confidence eiTors of the positions, of 
detection 1 and 2, respectively. Assuming that the probability 
distribution of a source position follows a Rayleigh distribution, 
the above formula (fill results in 99.73% completeness of the 
correlation (i.e. the probability, that we will detect all real 
counterparts, see e.g. [Watson et al.|2009 1. 
The catalogue was also compared to source catalogues in other 
wavelengths. In the near infrared we used the 2MASS catalogue 
(Skrutski e et al.|2006|l, while in the optical we used the LGGS 
catalogue ( Massey et al.|2006 l, for which we visually inspected 
the images by creating finding charts for all the X-ray sources as 
described in section|3.6| For the comparison with radio sources, 
we used the Braun (1990) catalogue (hereafter B90), based on 
Very Large Array (VLA) radio observations. The classification 
criteria are discussed in section |4] and summarized in Table |5] 
In order to find out which sources are associated with globular 
clusters in M31 we compared the catalogue positions with 
globular clusters from Galleti et al.| ( [2007 ) (and checked clas- 
sifications in the latest online version of the Revised Bologna 
Catalog (RBCv5, August 2012)). To identify X-ray sources with 
novae in M31, we correlated the catalogue with optical nova 
positions, obtained from the M 31 nova monitoring projecr] 



3.3. Variability 

To investigate the X-ray variability of the point sources, we de- 
veloped a Linux shell script, calling various CIAO tools, to ob- 
tain the source counts and upper limits for each source in each 
observation. We extracted source counts from the event files 
(with dmlist), using circular extraction regions. As radii, we 
used the 90% encircled count fraction PSF radius at 1 ke V, calcu- 
lated with mkpsfmap at the source position. Then we subtracted 
background counts calculated from the background maps ac- 
cordingly (with dmstat). Using the CIAO tool aprates we then 
calculated the source statistics. We executed two aprates runs 
with 68% and 99.7% confidence level for the errors and upper 
limits (coiTesponding to Icr and 3cr confidence for a Gaussian 
distribution). Luminosities were calculated from the aprates 



count rate output, as described in section 2. 1 



For each source, we calculated the fraction (ndet) of the num- 
ber of > Icr data points in the light curves over the number 
of observations, in which the source is in the detector's FOV. 
Additionally we introduced two measures of variability to the 
source catalogue (see Stiele et al.|2008[ Kong et al.|2002b| l. 
One is the variability significance (svar) 



SVar = (Fmax - Fir 



a)/yfc 



+ cr^ 



(2) 



where F is the flux and the cr are the errors associated with the 
flux. For the extreme values F^ax and F^in (fmax and fmin in the 
primary source catalogue table) we used the flux values, where 
F - cr is maximal and F + cr is minimal respectively (adopted 
from Sturm et. al. 2012, submitted to A&A). If a 3cr upper limit, 
lower than Fn,;,, was available, we used it as the minimum flux 
instead, and in place of its uncertainty cr„ 
used the Icr upper limit. 



in equation (J2]) we 



http://www.mpe.mpg.de/~m3 lnovae/opt/m3 1/index.php 
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Fig. 4. Population plot showing the source variability ratio fvar 
versus the maximum photon count rate of a source. Blue dia- 
monds are OBs, green squares are HVs, red circles are GlCs, 
black dots are novae and black crosses are the rest of the sources. 



The second variability parameter is the variability ratio (fvar), 
which stands for the ratio of maximum over minimum flux 
Fmax/Fmin of the sourcc. To calculate these properties we used 
data points for which the significance of the respective source 
given by aprates was above Icr. 

In addition to the overall variability significance svar between 
the maximum and minimum flux over 13 years, we calcu- 
lated the maximum of the variability significance (ssvar), be- 
tween neighbouring observations. The observations have differ- 
ent spacings in time but the maximum variability between them 
is still a good measure for variability on time scales from about 
ten days to one year. 

Using the data in Table 2, we plotted the HRC-I light curves for 
all the sources over the last thirteen years (visualized using mat- 
plotUb by Hunter 2007 1. The variability plots show brightness as 
count rates (label on right ordinate) and luminosity in M 3 1 (la- 
bel on left ordinate, assuming the power law spectrum from sec- 
tion 2. 1 and a distance of 780 kpc) with significance higher than 
Icr, as a function of time in modified Julian date (MJD). The 



as a function of time in modified Julian date (MJD). 
data points are filled circles (black) with Icr error-bars, show- 
ing the mean brightness during one observation. Whenever the 
source was outside the detector FOV of the relevant observation, 
there is no corresponding data point. If the source significance 
was below one, the 3cr upper limit of the source flux is marked 
by a red triangle. The observations span about 12.5 years from 
1999-1 1-30 to 2012-06-01. Sources No. 2, 4, 5, 8, 16, 22, 24, 44, 
52, 90, 303, 315 and 318 have no detection > Icr in the individ- 
ual observations and so we cannot show light curves for those 
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Fig. 5. Example of light curve plots showing X-ray luminosity 
(left axis, 0.2-10 keV) and photon count rate (right axis) over 
modified Julian date. Fluxes with Icr eiTors are given as black 
dots and 3cr upper limits as red triangles. 



There are 18 sources in the catalogue, for which the 90% encir- 
cled count fraction PSF overlaps with an adjacent source, thus 
causing some cross-talk (see "XT" flag in section [33| . As can 
be seen in Figs. [3] \T4\ [TS] and [T6] this overlap is mostly critical 
for the sources in the crowded core region of M 3 1 . 
All light-curve plots are available in the online version of this 
publication in Fig.|B.l| 
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3.4. Light-curve classes 

We derived two light-curve classes using the following selection 
criteria from the variability parameters defined in section[33| 



- Highly Variable (HV): ndet > 0.3 & (ssvar > 4 or svar > 5) 

- Outburst (OB): ndet < 0.3 & fvar > 15 & svar > 5 

HV sources are detected in > 70% of the observations and show 
strong variability between neighbouring observations and/or 
high variability at longer time scales in the HRC-I monitoring 
(e.g. No. 104inFig.|5]l. 

OB sources show a significant intensity increase by at least a 
factor of 15 and are not detected in at least 70% of the observa- 
tions, when in the field of view. We further divide OBs into three 
subclasses: 

- Short outbursts (SOB, e.g. No. 12 in Fig.|5]l, which decay to 
a quiescent level below the sensitivity limit in less than about 
500 days. 

- Long outbursts (LOB, e.g. No. 184 in Fig. [5]l, which need 
more than 500 days to decay. 

- Activity periods, which show strong variability, during an 
outburst in X-rays (AP, e.g. No. 185 in Fig.jSj. 

Some OBs may have been classified as transient X-ray sources 
in the literature, as they show one outburst and then disappear 
again. However, sources in M 3 1 are only detected if they are 
brighter than about 1x10^^ erg s~'. Therefore we cannot be sure, 
if a source is not still reasonably bright before and after OBs, 
sincelGarcia et al. 1(11998 ) find quiescent states of Galactic XRBs 

" • 1034, 



to be much lower (between 



10^^ erg S-' 



and 



erg s ) 



than the sensitivity limit. Good examples for such behaviour are 
No. 181 (Fig. |5]l or No. 152 and 198 in the catalogue, which 
show outbursts, but are still detected in between. Sources with 
very long APs can also be classified as HV instead of OB. 
Fig. |4] shows a scatter plot of the variability factor (fvar) versus 
the maximum photon count rate for all sources with variabil- 
ity information in the catalogue (created using TOPCAlrb. The 
plot shows a branching of variable sources at higher luminosi- 
ties. The OBs by definition have a higher fvar, the brighter they 
get, the HVs can have a low fvar at high luminosities (see also 
Fig. |7]). Due to the sensitivity limit of the observations, sources 
at lower luminosities generally have a smaller variability factor. 

3.5. Table description 

The catalogue consists of two tables, containing the 318 sources. 
The primary table with general source information (Table 1) and 
a secondary table with detailed source variability information 
(Table 2). 

The primary table (Table 1) contains the source number 
(SRCJD, Col. 1), the RA and Dec coordinates (J2000) in de- 
grees (RA, DEC, Col. 2, 3), three times the statistical Icr posi- 
tional error (from wavdetect, not including systematic errors 
due to centroid offsets far off-axi^ and the systematic cata- 
logue eiTor from offset correction and reference catalogue, see 
section |2.1| (eiTor, Col. 4), the net source counts and corre- 
sponding eiTor (NET_COUNTS, NET_COUNTS_ERR, Col. 5, 
6), background counts in the source region (BKG_COUNTS, 
Col. 7), the net count rate and corresponding error (NET_RATE, 
NET_RATE_ERR, Col. 8, 9), the background count rate and cor- 
responding error (BKG_RATE, BKG_RATE_ERR, Col. 10, 11), 

** http://www.starlink.ac.uk/topcat/ 

' http://cxc.harvard.edu/ciao/ahelp/wavdetect.html 



the exposure time of the most significant detection of the source 
(EXPTIME, Col. 12) and the highest source significance found 
in any observation (SRC_SIGNIFICANCE, Col. 13). Those pa- 
rameters are part of the wavdetect output. 
In addition we give the number of the highest resolution image 
the source is detected in (res. Col. 14, see section [TT] ), the en- 
ergy flux and its eiTor (flux, flux_err. Col. 15, 16), luminosity, as- 
suming the source lies in M 31 (see section 2.1 1 and error (lum, 
lum_err. Col. 17, 18) and the distance to the nearest neighbour- 
ing source (dist_NN, Col. 19). 

From the light curves of each source we give the parameters svar, 
fvar, ndet, ssvar, fmax, fmin and errors (Col. 20 to 27) and the 
derived light curve classification (lc_class. Col. 28), described in 
When a source is detected only in the merged im- 



3.3 



section 

ages or only once in the detector FOV, we cannot give any light 
curve information and the variability values are set to 0.0. 
We list the X-ray correlations of the sources with other X-ray 
catalogues, described in 3.2 (XID, Col. 29). Due to the high 
X-ray source density in the central area of M31, there can be 
confusion of X-ray sources, especially in the case of catalogues 
obtained with lower resolution instruments, since they do not re- 
solve all sources. The column contains a string of correlations, 
which provides the source No. in the catalogue, followed by the 
name of the reference catalogue, the source number or name and 
if available the type of the source, correlating with the catalogue 
source (e.g. [132, 'PFH2005', ' 310 XRB; ']). 
Names of optical novae and globular clusters, associated with an 
X-ray source are listed separately (nova, glc. Col. 30, 31). For the 
globular clusters we also list the quality flag of the RBCv5 (1: 
confirmed, 3: controversial object, 6: confirmed star, see GaUeti 
iCt al. 2007 ), stating if the globular cluster has been confirmed. 
Source No. 85 in the catalogue correlates with a former GIC 
candidate and has since been confirmed to be a foreground star. 
Source No. 49 correlates with a controversial object from the 
RBCv5 and we classified it as GIC candidate. 
We then give the source type (class. Col. 32), based on the iden- 
tifications in the X-ray, radio and optical wavelengths, as de- 
scribed in section 3.2 and the criteria summed up in Table l5] 
Candidates are designated by <>. 

If necessary, we comment (comment. Col. 33) on special proper- 
ties of a source. Comments refer to special light curve properties 



of a source (see section 3.4 1 or note the lack of a light curve (no 
Ic) for some sources which are only found in the merged images, 
or too far off the optical axis (we did not include data points, 
where the extraction region overlaps with the detector edge, see 
For some sources we added a bad position flag 



3.31 



section 

(bad position) in the comments, since we found that the CIAO 
tool wavdetect did not find the centroid position for some faint 
and far off-axis sources. If the identification with known X-ray 
sources or counterparts in radio or optical data, due to the bad 
position is not clear, we do not provide a classification for the 
sources based on these identifications (in Col. 32). If the 90% 
encircled count fraction extraction region of a source overlaps 
with another extraction region this is highlighted by "XT" fol- 
lowed by the overlapping source No. (e.g. XT[164]). For each 
statement in the "comments" column, which is not derived from 
our analysis, we give an abbreviation of the relevant reference. 
The time variability table (Table 2) contains the observation 
numbers (ObsID, Col. 1), start date of the observation (date, 2), 
the modified Julian date at the start of the observation (MJD, Col. 
3) and three columns for each source, giving its flux and upper 
and lower 68% confidence limits (Col. 4 to 957). If a source is 
not detected > Icr, we give a 3cr upper limit in the upper confi- 
dence limit column and zero for flux and lower confidence limit. 
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Table 4. Prominent optical identifications up to about 1" separation from one of the X-ray sources. 



No." 


RA(J2000) 


Dec(J2000) 


XID* 


OptID ' 


V 


B-V 


U-B 


V-R 


R-I 












[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


26 


00:42:10.30 


+41:15:10.2 




J004210.30+411510.6 


19.31 


1.80 


1.20 


0.69 


1.12 


37 


00:42:15.54 


+41:20:32.0 




J004215.52+412031.5 


20.52 


0.42 


-0.47 


0.43 


0.83 


60 


00:42:26.18 


+41:25:51.1 




J004226. 17+412552.0 


20.61 


1.09 


-0.80 


0.23 


0.44 


75 


00:42:32.07 


+41:13:14.6 


<XRB> 


J004232.08+411315.2 


20.75 


0.92 


0.74 


0.48 


0.57 


85 


00:42:33.48 


+41:21:38.9 


fgStar 


J004233.56+412138.5 


17.18 


1.61 


1.13 


1.11 




94 


00:42:35.18 


+41:14:21.4 




EO 












185 


00:42:45.93 


+41:10:36.9 


<XRB> 


J004245.84+41 1036.9 


21.37 


0.43 


-0.73 


0.41 


0.51 


207 


00:42:49.06 


+41:19:46.1 




J004249.07+41 1945.8 


19.81 


0.63 


-0.02 


0.84 


0.60 


238 


00:42:55.18 


+41:18:36.3 


<XRB> 


EO 












262 


00:43:01.74 


+41:17:27.3 


< fgStar > 


J004301. 78+41 1726.5 


18.05 


1.16 


0.80 


0.72 


0.75 


264 


00:43:03.03 


+41:20:41.2 




EO 












283 


00:43:13.16 


+41:18:13.3 




J004313.18+411814.1 


22.19 


1.35 




0.26 


1.24 


284 


00:43:14.34 


+41:16:50.8 




EO 












289 


00:43:17.00 


+41:12:24.9 


< fgStar > 


J004317.04+41 1224.2 


18.75 


1.59 


1.02 


1.00 


1.28 


291 


00:43:17.83 


+41:11:13.4 


< fgStar > 


J004317.91+411113.5 


19.04 


1.67 


1.10 


1.08 


1.48 


292 


00:43:18.62 


+41:09:49.8 


nova 


J004318.58+410949.1 


21.34 


0.18 


0.01 


0.41 


-0.23 


294 


00:43:18.91 


+41:20:16.7 


<SSS> 


EO,J004318.89+412017.2 


22.43 


0.70 


-1.45 


0.49 


0.94 


301 


00:43:24.81 


+41:17:26.5 




J004324.81+41 1726.1 


21.74 


-0.00 


-0.79 


0.02 


-0.72 



Notes. *"' X-ray source number in the catalogue. '*' X-ray class of the source in the catalogue. ''' Source name of optical identification with a 
source from.Massey et al.,^2006j, or EO for emission line objects, found while inspecting the images. 



If a source is out of the FOV the three flux columns are zero. 
Sources, which are not included in previous dedicated M 31 cat- 



alogues are listed in Table 10 For the new sources as for all other 
sources in the catalogue, we give a light curve classification (see 
section 3.4 1 and hardness information if possible. For some faint 



sources we can only give the positions since we do not have any 
significant light curve information. 

3.6. Optical counterparts 

The search for optical counterparts is important to further con- 
strain the nature of an X-ray source. We checked finding charts 
derived from broad and narrow band images, provided by the 
Local Group Galaxy Survey project (LGGS, see Masseyetal. 
2006[l. Where we found a prominent and not yet clearly identi- 



fied optical source within the 3cr error circle and up to a separa- 
tion of about 1" of the catalogue source, we list the source name 
from the jlVrassey et al.| (20Q6J catalogue, with visual magnitude 
and optical colours in Table l4] We regard an optical counterpart 
as prominent, if we visually identified it as the brightest and at 
the same time closest counterpart in the LGGS images. Table 
ffllists the closest LGGS counterpart to the corresponding X-ray 
position. Expected counterparts in this list could be galactic fore- 
ground stars, background galaxies or high-mass stars (as com- 
panion stars in high-mass X-ray binaries in M 31, see |4.1.5| l. We 
do not include in this list confirmed objects from earlier studies 
such as galaxies, foreground stars, globular clusters or supernova 
remnants. 



4. Discussion 

In the central field of M31 there are many X-ray sources, 
which have previously been classified or identified into different 
classes according to their X-ray spectra or at least HRs, their 
time variability and optical and radio correlations (see X-ray 
catalogues in 3.1 1. Since HRC-I has only limited spectral reso- 
lution, we had to rely on previous catalogues or neighbouring 



Chandra ACIS, Swift XRT or XMM-Newton EPIC observations, 
to estimate the hardness of a source (see section [23] l. 
New in this catalogue, however, are the good source positions 
in the inner area of the field due to the high resolution of the 
HRC-I detector. The long integrated exposure improves the 
detection of faint sources. Table [5] gives an overview of the 
criteria for different classifications in the catalogue, which are 
discussed in more detail throughout this section. Table l5] also 
lists the statistics for how many sources are clearly identified in 
the catalogue and how many are classified using the selection 
criteria. 

For most sources we were able to analyse time variability in 
detail. This enabled us to give a classification based on light 
curves, which has not been possible to this extent in previous 
catalogues. However, there are still sources without classifica- 
tion in the catalogue, because we do not have significant light 
curve information on them and thus can not obtain a light curve 
class (see section [33| ). 

In general, the variability, which we infer from the light curves 
can be caused by intensity or spectral changes, resulting from 
anisotropic radiation of a rotating object, full or partial occul- 
tation of the X-ray source or a change in energy output of the 
emitting source. 

In the catalogue we classify 195 sources as located in M31 
(61.3 %) and 13 foreground and background sources (4.1 %). 
110 sources remain unclassified (34.6 %, see Table [5]l. Details 
of the classifications are given below. We regard a source as 
identified if the source has been clearly identified in previous 
studies and if this identification is compatible with the HRC-I 
analysis. Otherwise we regard a source as classified, if it has 
only been a candidate in previous studies or if we get the 
classification only from the variability analysis (candidates 
indicated by <> in the catalogue). 
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Table 5. Source classes and criteria in the catalogue. 



Class " 



Criteria 



Identified Classified 



XRB hard X-ray source, HV, OB 

GIC correlation with optically confirmed GIC, hard spectrum 

SSS X-rays mostly below 1 keV 

nova correlation with optical nova 

SNR softer X-ray source and/or extent, correlation with radio or optical SNR 

AGN hard X-rays, radio identification 

fgStar softer X-rays, optical identification 

unclassified 



3 


101 


32 


1 





6 


45 




6 


1 


1 


1 


7 


4 


110 





Stiele et al. 



Pietsch et al. 



Notes. *"* Classes are adopted from 

nova: optical nova, SNR: supernova remnant, AGN: active galactic nucleus, fgStar: foreground star. 



(2011 1; 



1 2005b^. XRB: X-ray binary, GIC: globular cluster, SSS: super-soft source. 



4.1. X-ray binaries (XRB) 

According to theoretical models. X-ray binaries (XRB) are 
thought to consist of a compact object (black hole, neutron star 
or white dwarf) and a donor star of various types, from which 
matter is accreted onto the compact object. 
From the phenomenological point of view, we divide X-ray bi- 
naries into globular cluster (GIC) sources (optically identified), 
XRBs in the field with a hard X-ray spectrum (probably a black 
hole or a neutron star as the compact object) and XRBs with a 
super-soft (SSS, nova) X-ray spectrum (probably a white dwarf 
as the compact object). 

While for hard sources a significant part of the photons are de- 
tected at energies > 1 keV, for SSS sources the majority of pho- 
tons are detected below 1 keV. 



4.1.1. Globular cluster sources (GIC) 

We find 32 of the X-ray sources in the catalogue to be asso- 
ciated with confirmed M31 GlCs (excluding novae). From the 
variability criteria (see section 3.4 1, we classify 16 of them as 
highly variable (HV, e.g. No. 72, shown in Fig.l6]l and 2 as OB. 
14 GIC sources have no variability class, because they are too 
faint in X-rays to detect significant variability. In addition there 
are two GIC sources associated with optical novae (No. 81 and 
150 in section [4. 1.6[ ). The eiTors of both the X-ray positions and 
the globular cluster catalogue (section [3T| are small and the den- 
sity of globular clusters in the central region of M 3 1 is too low 
to lead to a significant contribution of chance coincidences. This 
allows us to state with high confidence, that the X-ray sources, 
we find correlated with GlCs are actually inside an M 31 globu- 
lar cluster. 



le37 source 72 (00:42:31. 26,+41:19:38.9) ie-2 



Clark] ( |1975| l and [Fabian et al] ( |1975| l found, that bright X-ray 
sources in GlCs in the Milky Way are LMXBs and that the den- 
sity of LMXBs per stellar mass in GlCs is higher, than in the 
field. They explain this fact by two different dynamical mecha- 
nisms of binary formation (three body scattering and tidal cap- 
ture). Based on these findings, we assume, that all the GIC X-ray 
sources are LMXB candidates. Since all the GIC sources with 
significant detections in their light curves show high variabil- 
ity between neighbouring observations and are thus classified as 
HV or OB sources, we classify all sources in the field showing 
the variability class HV or OB as XRB candidates. As shown in 
Table|5]and the following sections, other source types can either 
be excluded due to time variability behaviour, spectral informa- 
tion or identifications at other wavelengths. 
The classification scheme gets further support by the fact, that 
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Fig. 6. Light curve plot of the GIC X-ray source No. 72, showing 
X-ray luminosity (left axis, 0.2-lOkeV) and photon count rate 
(right axis) over modified Julian date. 



similar time variability in the 0.2 to 10 keV X-ray energies has 
been observed in Galactic X-ray sources, which have been con- 
firmed to be XRBs (e.g. Homan et al.r2010). The variability in 
the HRC-I light curves might be enhanced by spectral changes 
of a source. 

As mentioned above, the GIC sources which have not been 
classified by the light curve variability are most probably also 
LMXBs. Therefore, also in the field, several LMXBs will not be 
classified by the light curve criterion. 

There are no other source types, besides the SSS emission of op- 
tical novae in M31, that correlate with the OB and HV sources 
and their distribution becomes more dense toward the centre of 
M31, following the density of stellar populations in the galaxy. 
Thus we are very confident, that the OB and HV sources are 
actually XRB candidates located in M 31, as we would expect 
background or foreground objects to be randomly distributed 
across the FOV. 



4.1.2. XRB variability 

Fig. It] and Fig. ISlshow the spread of X-ray luminosity Lx, we 
observe for HV and OB XRB (HVXB, OBXB) sources in the 
field of M31. The OBXBs all show a high variability factor 
(as seen in Fig. Hli. Several sources are detected in a very deep 
observation at a flux below the most stringent upper limits. 
Most HVXBs do not fall below the sensitivity limit during the 
observed period of thirteen years, but show a more confined Lx 
range. Some of them have a particularly low Lx spread, which 
could hint at an underlying source class within the HVXBs. 

Phenomenological discussion of some XRB candidates, 
with hard X-ray spectrum: 
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Fig. 7. Variability plot showing the spread of luminosity (left y-axis) and count rate (right y-axis) for all HV XRB sources (HVXB) 
in the field of M 31. The x-axis is labelled with the catalogue source number (No.). If the minimum count rate is derived from an 
upper limit of the light curve, the lower count rate is marked by a downward arrow. 
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Fig. 8. Variability plot showing the spread of lu- 
minosity (left y-axis) and count rate (right y- 
axis) for all OB XRB sources (OBXB) in the 
field of M31. The x-axis is labelled with the 
catalogue source number (No.). If the minimum 
count rate is derived from an upper limit of the 
light curve, the lower count rate is marked by a 
downward arrow. 



No. 3 shows an AP, starting in 2004 (see Trudolyubov 
et al.|2006| l. It was first detected by |Galache et al.| ( 2007af 



No. 47 shows a SOB starting on 2009-11-07, followed by an 
AP lasting until the end of the HRC-I monitoring in 2012. In 
an ACIS-I observation (ObsID 11275 from 2009-11-11), the 

source shows a hard spectrum (see Table [3)l. 

No. 49 is an X-ray burster (see Pietsch & Haberl 2005 1, showing 
a highly variable light curve. The source is associated with a 
GIC candidate (see RBCv5 in section [J!2| ). 
No. 54 shows two high luminosity states, each lasting for several 
years. During the low state (luminosity around 5x10-''' erg s"'), 
the variability is low and during the high state (luminosity 
around 2x10^^ erg s"'), variability is high. From 2008 until 
2010-02-15 (ObsID 9825) the source steadily brightens until it 
reaches the high state. During this high state, the source shows 
variability similar to AP sources. This system is discussed as 
being a black-hole X-ray binary, showing spectral state changes 

in |Barnard et~aL] ( |20n] ). 

No. 57 shows two APs in 2004 ( |Stiele et al.| '2008') and 2007 
( [Henze et al.|2008a) . The source is not bright enough to classify 
it as OB (see criteria in section 3.4 1. We still classify the source 



as XRB candidate, based on the literature. 

No. 58 is not confirmed as a fgStar candidate as in Stiele et al. 

( |2011| l, as there is no star visible in the Massey et al. (2006i 

images and we find a highly variable light curve. Therefore we 

classify the source as an XRB candidate. 

No. 65 shows a SOB, lasting a few months and starting on 2008- 

09-01 (in ACIS-I ObsID 9523) with a hard X-ray spectrum (see 

Table|5]l. In an observation on 2008-07-15 (ObsID 9522) it was 

not detected. It is further detected with ACIS-I on 2008-10-13 

(ObsID 9524) and not detected after 2008-11-27 (ObsID 9521), 

which constrains the duration of the outburst. 

No. 67 shows an AP s tarting 2009-05-29 (see Galach e et al.| 

2009HKong et al.|2009| . The source was active untill 2011-03- 



10 in HRC (ObsID 13180). 

No. 82 shows a SOB in the globular cluster B112, starting 

2012-06-01 (Hofmann et al. 2012;). In following monitoring 



observations with the Swift XRT, the source showed a rapid 
decay within a few weeks. It is one out of only two OB sources 
we found in globular clusters (see section [4. 1.2 1. 
No. 84 is a very faint source in the catalogue. Williams et al. 
(2005c) detect the source in an outburst in their ACIS-I monitor- 
ing program and discuss it as a recuiTent transient, since it was 
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already detected by the ROSAT HRI instrument (see |Primini| 

|etal.|1993| l. 

No. 104 is the most luminous HV source in the catalogue. It has 

been discussed as a Z-track XRB by |Barnard et al. (2003). 

No. 116 shows a SOB, which is only detected with HRC-I, 

but not discovered with ACIS. We extract 34 counts in the soft 

HRC-I band (PI channel 50 to 130), 69 counts in the medium 

band (131 to 200) and 25 counts in the hard band (201 to 280). 

This count distribution hints at a hard X-ray spectrum (see 

IPietsch et al.|2007| . 

No. 127 shows a SOB in an HRC-I observation starting 2006- 

11-13. ACIS-I ObsID 07064 (2006-12-04) indicates a hai'd 

X-ray spectrum (see Table[3]l. 

No. 128 shows an AP lasting from 2006-12-04 (ACIS-I ObsID 

7064), until 2010-02-04 (ACIS-I ObsID 11278) with a hard 

X-ray spectrum (see Table[3]l. 

No. 132 shows a LOB, which we cannot clearly constrain, 

since the source is already bright in the first HRC-I observation. 

It was discovered by Murray et al. ( 1999| l in October 1999. 

Trudolyubov et al. (2001|l argued, that the source shows an 



outburst of at least 400 days, with a plateau phase of about one 
year, which should indicate a black-hole binary. 
No. 149 is the second ultra-luminous X-ray (ULX) source, 
detected in M31 (see |Henze et al. 2012b a). It was also ob- 



2013 1 and in radio 



served with Swift XRT (see Esposito et al. 

observations by|Middleton et al.|(]20131l. 

No. 153 and 158 are clearly resolved as two sources in the 

catalogue (see also Kong et al.|2002b| ). In |Barnard et al. ( 2012a| l 

the source is treated as one source. We give two light curves, 

which show some cross talk. 

No. 161 shows a SOB and is first detected in ACIS-I ObsID 

12973, starting 201 1-08-25. In the following ACIS-I observation 

(ObsID 12974) from 2011-09-28 the source can be classified as 

hard (see Table [3]). 

No. 181 is reported as an X-ray transient by Shirey (2001 ). In 

the HRC-I monitoring light curve, with lower sensitivity limit, 

we detect the source as HV. 

No. 185 shows an AP from 2008-05-31 (ACIS-I ObsID 9529) 

till 2010-03-05 (ACIS-I ObsID 11279). In HRC-I observations 

we detect the source as HV. The source has a hard spectrum, 

derived fmm XMM-Newton observations ( |Pietsch et al.|2008bl l. 

No. 201 shows a long-term variability, with two LOBs at the 

beginning and end of HRC-I observations. After the second 

outburst, the source stays in a high-luminosity state, which 

slowly declines over the course of three years and is still at 

3x10^^ erg s ' at the end of the HRC-I monitoring. 

No. 222 shows a SOB in just one of the HRC-I observations 

(ObsID 7284 on 2006-09-30). ACIS-I ObsID 7064 on 2006-12- 

04 indicates a hard spectrum (see Table l3]l 



No. 234 is the first ULX source found in M 31 (see Henze et al. 



2009a 



2012a 



IKaur et al.poTIl [Middleton et al.|[20T2l [Barnard et al 



|Nooraeeetal.|2012| i. 



No. 237 is a HV source. In Nov. 2010 it dropped in luminosity 
from ~ 10^** erg s"' and after that showed high variability. This 
indicates a spectral state change. 

No. 240 shows rather slow variability. It is not classified as HV 
as the count rate has large errors, when the source is brightest. 
No. 247 has a SOB starting 2010-10-19 ( Barnard et al.||2012a| l 
in ACIS-I ObsID 12160, where we find a hard spectrum (see 
Table[3]l. 

No. 259 shows a faint AP from 2006-06-05, until 2010-12-12, 
which is not significant enough for the outburst criteria. ACIS-I 
observation ObsID 7064 on 2006-12-04, indicates a hard 
spectrum (see Table l3]l. 



No. 275 is an X-ray burster in a globular cluster system, which 

is discussed by Pietsch & Haberl (2005 i. 

No. 276 shows dips with a 107 min X-ray period (Mangano 



|et al.| 2004| l and was thus classified as an X-ray binary. In 
the HRC-I light curve the source does not show high enough 
variability to be classified as HV. 

No. 278 is at the position of a confirmed globular cluster at the 
edge of a SNR in the optical images (see Massey et al.||2006 1, 
which might cause some contamination of the light curve of the 
GIC source. However, the SNR has to be very faint in X-rays, 
since we do not see an extended source at the position in the 
merged X-ray images 



No. 279 shows a AP starting 2004-05-23, discussed in Williams 

|et al. (2005b i. The last detection in the HRC-I light curve is on 

2005-02-21. 

No. 285 shows a very peculiar light curve with two rather well 

separated luminosity states and has previously been found to 



have a 2.78 h period by |Trudolyubov et aLlP002| l. [Barnard et al. 
( 2006 1 argue, that there is another component to the variability, 
caused by precession of the accretion disc of the compact object. 
No. 296 shows a SOB, starting 2008-08-21, discovered with 
Swift XRT (Pietsc h"eni.||2008a[palache et al.||2008c| [Heiize 
,et al.,2008b,,Galache et al.|2008a| l. 

No. 299 is a HV source. One of its brightness peaks has been 
reported by'Nooraee et al. ( 2009 1. 

No. 300 shows one SOB and is only detected with HRC-I. To 
estimate the hardness of the spectrum, we extract 629 counts 
in the soft HRC-I band (PI channel 50 to 130), 1440 counts in 
the medium band (131 to 200) and 479 counts in the hard band 
(201 to 280). This distribution indicates a hard X-ray spectrum 
dPietsch et al.|2007| ). 

Recurrent outburst hard X-ray sources: 

No. 15 shows three APs in the light curve, the first of 
which is discussed in Osborne et al.||2001) and [Williams et al 



(2006|. The source is discussed as a recurrent transient in Stiele 



|et al. ( |2011| l, since the source was already observed with the 
ROSAT observatory (see [Primini et al. 1993| l. We detect the 
latest AP, with the source still active at the end of the HRC-I 
monitoring on 2012-06-01. 

No. 40 shows many recurrent SOBs about every 1 .2 years (see 
[Henze et al.|2007) [Pietsch et al.|2010at|Galache et al.|2007b|l. 
No. Ill shows recurrent SOBs. One is discussed in lGarcia et al.l 
( |200T] l. The source was not active in ACIS-I 9521 (2008-11-27) 
and shows one outburst in HRC-I ObsID 10684 (2008-12-18). 
We detected a second SOB starting in ACIS-I ObsID 12164 
(2011-02-16), which we also observe one day later in HRC-I 
ObsID 13178. If there is a quasi-periodic recurrence time, the 
HRC-I data would suggest an outburst period of ~2 years. 
No. 140 shows recurrent SOBs and a hard spectrum in ACIS-I 
([ Kong et aL][2002b| l and XMM-Newton EPIC observations. 



separated by ~6 years. 
No. 147 shows either short recurrent OBs about every year, 
which are missed in some of the HRC-I monitoring observations 
or it shows two APs. 

No. 162 shows a SOB in HRC-I ObsID 7283 (2006-06-05). 
ACIS-I obser-vation (ObsID 7138 from 2006-06-09) closely 
after the HRC-I detection (ObsID 7283 from 2006-06-05) 
indicates a hard spectrum (see Table |3]l. The earliest detection 
of the outburst was visible in an ACIS-I observation from 
2006-05-26 with a total of 9 counts. In earlier observations there 
was an outburst only seen with ACIS-I in an observation starting 
on 1999-10-13 (ObsID 303), where we extract 26 counts. Due 
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Table 6. Additional Chandra detections for recurrent outburst Table 7. Additional detections for recuiTent outburst source No. 
source No. 246. 290. 



Instrument 


ObsID" 


Date* 


ACIS-I 


303 


1999-10-13 


ACIS-I 


312 


2000-08-27 


ACIS-S 


1581 


2000-12-13 


ACIS-S 


1854 


2001-01-13 


ACIS-S 


1575 


2001-10-05 


ACIS-I 


7136 


2006-01-06 


ACIS-I 


7064 


2006-12-04 


ACIS-I 


8183 


2007-01-14 


ACIS-I 


8191 


2007-06-18 


ACIS-I 


8192 


2007-07-05 


ACIS-I 


9523 


2008-09-01 


ACIS-I 


10715 


2009-09-18 


ACIS-I 


10716 


2009-09-25 


ACIS-I 


11838 


2010-05-27 


ACIS-I 


12161 


2010-11-16 


ACIS-I 


12162 


2010-12-12 


ACIS-I 


12973 


2011-08-25 


ACIS-S 


14197 


2011-09-01 


ACIS-S 


14198 


2011-09-06 



Notes. *"* Observation identification number. **' Start date of the obser- 
vation. 



to this outburst the source is part of the |Kong et"ar ( 2002b| l 
catalogue and thus is classified as a source with recurrent 
outbursts. 

No. 189 shows recuiTent SOBs in the HRC-I light curve. The 
source brightens from 2009-11-07 till 2010-02-26. ACIS-I 
observation ObsID 11279 on 2010-03-05 indicates a hard 
spectrum (see Table [3]l. The source is active till ACIS-I ObsID 
11839 from 2010-06-23. 

No. 246 shows multiple SOBs, for which we catch three in a 
very bright state in the HRC-I observations. There are many 
additional observations with the Chandra ACIS detector (see 
Table l6|l, filling in some of the HRC-I observation gaps. The 
source seems to have bright outbursts about every year, lasting 
up to about three months. 

No. 290 shows OBs, lasting up to several months, about every 
year during the time span covered by HRC-I observations. Table 
[TJlists additional observations from other instruments, where the 
source is visible in the X-ray images, which fill in some of the 
time, covered by HRC-I hght curves. The XMM-Newton EPIC 
HR of the source indicates a hard X-ray spectrum (see 2XMM 
catalogue by Watson et al.|2009]l. 

No. 313 shows recurrent OBs, first discovered with the Einstein 
and ROSAT X-ray observatories. In XMM-Newton EPIC 
observations it showed a harder X-ray spectrum (see 2XMM 
catalogue by [Watson et"aL]|2009| l. In the HRC-I light curves 
we detect four OBs. It shows both LOBs and SOBs, with the 
SOBs being about nine months apart. Table l8] gives additional 
observations, where the source is detected. 



4.1.3. GIG population 

We found that after applying the variability criteria for XRB 
selection, the ratio of hard XRBs classified as HV (HVXB) to 
hard XRBs classified as OB (OBXB) is 2.3/1 in the field of 
M 31. For the hard XRBs found in GlCs we obtained a ratio of 



Instrument 


ObsID" 


Date* 


ACIS-S 


1582 


2001-02-18 


ACIS-S 


2900 


2002-11-29 


XMM-EPIC 


0405320701 


2006-12-31 


XMM-EPIC 


0405320801 


2007-01-16 


XMM-EPIC 


0405320901 


2007-02-05 


ACIS-I 


8184 


2007-02-14 


ACIS-I 


8185 


2007-03-10 


XMM-EPIC 


0505720301 


2008-01-08 


XMM-EPIC 


0505720401 


2008-01-18 


XMM-EPIC 


0505720501 


2008-01-27 


XMM-EPIC 


0505720601 


2008-02-07 


ACIS-I 


10553 


2009-03-11 



Notes. *"' Observation identification number. **' Start date of the obser- 
vation. 



Table 8. Additional observations for recurrent outburst source 

No. 313. 



Instrument 


ObsID" 


Date* 


ACIS-S 


309 


2000-06-01 


ACIS-S 


1581 


2000-12-13 


ACIS-I 


4680 


2003-12-27 


ACIS-I 


4681 


2004-01-31 


XMM-EPIC 


0405320501 


2006-07-02 


XMM-EPIC 


0405320601 


2006-08-09 


XMM-EPIC 


0405320701 


2006-12-31 


ACIS-I 


8183 


2007-01-14 


XMM-EPIC 


0405320801 


2007-01-16 


XMM-EPIC 


0405320901 


2007-02-05 


ACIS-I 


8184 


2007-02-14 


ACIS-I 


8185 


2007-03-10 


ACIS-I 


7068 


2007-06-02 


XMM-EPIC 


0505720501 


2008-01-27 


XMM-EPIC 


0505720601 


2008-02-07 


XMM-EPIC 


0650560201 


2010-12-26 


XMM-EPIC 


0650560301 


2011-01-04 


XMM-EPIC 


0650560501 


2011-01-25 


XMM-EPIC 


0650560601 


2011-02-03 


ACIS-I 


12164 


2011-02-16 


ACIS-S 


14197 


2011-09-01 


ACIS-S 


14198 


2011-09-06 


XMM-EPIC 


0674210201 


2011-12-28 


XMM-EPIC 


0674210301 


2012-01-07 


XMM-EPIC 


0674210401 


2012-01-15 


XMM-EPIC 


0674210501 


2012-01-21 


XMM-EPIC 


0674210601 


2012-01-31 



Notes. '"' Observation identification number. '*' Start date of the obser- 
vation. 



7.5/1 (HVXB/OBXB). For this statistical analysis we only used 
sources inside the ~ 70% effective exposure radius (see Fig. |2]i, 
to avoid effects by the inhomogeneous coverage outside this cir- 
cle. Fig.l9]visualises, that most GIC sources are HVXB sources 
with rather low variability factor. This can also be seen in Fig. 
|4] where the GIC sources concentrate in the HV source branch. 
Since we have fewer sources (and thus higher statistical uncer- 
tainty) in the GlCs, we estimated the probability that this dif- 
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Fig. 9. Variability plot showing the spread of luminosity (left y- 
axis) and count rate (right y-axis) for all OB or HV GIC sources. 
The X-axis shows the source No. in the catalogue. If the mini- 
mum count rate is derived from an upper limit of the light curve, 
the lower count rate is marked by a downward arrow. The blue 
diamonds mark the sources which are classified as OB sources 
in the catalogue. 



ference could occur by chance. We used Binomial statistics and 
assumed that the true ratio is 2.3/1 (from the field sources, where 
we have a sample of 97 sources, 29 OBXB and 68 HVXB). The 
probability to obtain a ratio of 7.5/1 after 17 trials (2 OB and 
15 HV sources in GlCs) is only 8%. This may indicate, that the 
abundance fraction of outburst sources in GlCs might be lower 
than in the field. The reason for this possible difference could be 
differing formation processes of LMXBs in the GIC and field en- 
vironments. However, there is also the possibility of additional 
OBs in GlCs, masked by interfering HV sources, which could 
not be resolved if they were located in the same GIC. A higher 
OB fraction in the field might also be caused by high-mass X-ray 
binaries (HMXBs, e.g. Be-XRBs with Type II outburst, see |Reig| 



201 l] l in the sample, which were not optically identified. It has 
be en discussed by|Clark|( I975| l, Fabian et al. ( 1975 1 and recently 



by |Voss & Gilfanovl ( |2007bf that XRBs in GlCs are likely to 
form through tidal capture in the dense star population of a glob- 
ular cluster, which they also use to explain an over-abundance of 
XRBs in GlCs. 

Overall, the ~ 70% effective exposure radius of the catalogue 
contains 154 confirmed GlCs from the RBCv5 (see section [372] i. 
These are associated with 27 X-ray sources, suggesting that at 
least 17% of GlCs in M31 should host an XRB. This is a very 
conservative lower limit, since the observations of M 3 1 are not 
deep enough to detect all XRB sources. 

4.1 .4. Low mass X-ray binaries (LMXB) 

Following the variability analysis, the majority of X-ray sources 
in M 31 seem to be LMXBs (see section [4. 1.2[ ). If one of those 
variable sources shows a hard X-ray spectrum, we classify it as 
an XRB candidate (< XRB >). The HRC-I observations mainly 
cover the bulge of M 31, which hosts a predominantly old stellar 
population. Therefore, we would only expect very few HMXB 
candidates in the catalogue (for an estimate see |Voss & Gilfanov] 
[2007a). 

Whether the compact object is a neutron star or a black hole can 
to some extent be decided from X-ray spectra or some variability 



properties (soft X-ray sources with transient outbursts are most 
likely black hole systems, see TanakapOOO i. The variability of 
XRBs in the light curves could in part be due to the transition 
from one spectral state into another, where hardness ratios in the 
0.2-10 keV X-ray band would change. Such state changes in X- 



ray binaries have been discussed by Hasinger & van der Klis 
(1989). For similar results from the RXTE satellite see |Homan 
et al . (2010|. 

The Eddington luminosity (proportional to the mass of the com- 
pact object in XRBs) describes the limit at which the outward 
photon pressure becomes so high, that no further material can be 
gravitationally accreted. So from the luminosities at which dif- 
ferent spectral states of a source exist, a rough estimate for the 
mass of the compact object can be derived (see discussion of No. 
54 in |Bamard et al.|201 l| l. If the mass is between about 1.4 solar 
masses (Mq) and 3.0 Mq, we would expect a neutron star. If the 
mass exceeds 3.0 M© the compact object should be a black hole 
(see |Tanaka|2000) . 

Our sample contains 4 sources exceeding the Eddington limit for 
a 3.0 Mq object. Two of those are the ULX sources, which have 
already been discussed as BH XRBs in previous studies (e.g. 
|Kaur et a l. 2012; Middleton et al. 2013). From HRC-I data alone 
a clear distinction between BH and NS systems was not possible 
since there was not sufficient simultaneous spectral information. 
Short bursts (at time scales of ~ one minute) or pulsations in the 
short-term light curves (during an observation) of a source are 
strong evidence for a neutron star binary. Bursts (as have been 
found for No. 49 and No. 275 by [Pietsch & Haberll[200 5l) are 
thought to arise from short nuclear burning flashes on the surface 
of a neutron star. For a recent review on neutron star binaries see 



Bhattacharyya]([20T0) l. 

Garcia et al.| ( 1998 1 discuss the separation of some galactic neu- 
tron star and black hole binaries, using the variability range of 
OBs. The sensitivity limit of the HRC-I observations for sources 
in M 3 1 does not allow to measure high enough ratios between 
maximum and minimum X-ray luminosity to apply the criteria, 
proposed by Garcia et al. ( I998|l, to the sample of sources. 



4.1.5. Higli mass X-ray binaries (HIVIXB) 



HMXBs are expected in star forming regions, which are mainly 
located in the disk of M3I. So far there are no confirmed 
HMXBs in M31. We would classify a source as HMXB 
candidate if it shows XRB-like time variability in the light 



cui-ves (see criteria in section 4.1.2i, has a hard X-ray spectrum 



and shows an optical counterpart with colours consistent with a 
high-mass star in the Massey et al. (2006 1 images. 
Source No. 301 is bright in X-rays, has a prominent optical 
counterpart (LGGS and Hubble space telescope (HST) images) 
and shows no radio emission in the reference catalogues. The 
colours and magnitude of the optical counterpart (see Table 
|4| and the hardness ratios of the X-ray emission would be 
consistent with a HMXB in M31 (see Sturm et al. 2012, 
submitted to A&A, for a comparison with HMXBs in the Small 
Magellanic Cloud). However the errors of the X-ray hardness 
and the optical colours are too large to exclude an AGN as 
the source of the X-rays. We do not detect Ha emission at the 
source position and have no spectrum of the optical counterpart. 
The HRC-I X-ray light curve is too faint to show significant 
variability. 
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Fig. 10. Light curve plots of LOB X-ray sources associated with 
optical novae, showing X-ray luminosity (left axis, 0.2-lOkeV) 
and photon count rate (right axis) over modified Julian date. 



4.1.6. Super Soft Sources (SSS) and Novae (nova) 

The class of SSS in M31 is thought to be associated with 
cataclysmic variables (CV), which are binary systems with a 
white dwarf as the accreting object. Those systems are also 
responsible for nova explosions (see jPietsch et al.|2007| . Most 
of the sources with super-soft hardness ratios (photon counts 
mostly below 1 keV) are associated with nova explosions, which 
are discussed in detail in Pietsch et al.| (p005a|, |Pietsch et al.| 
(2007), Henze et al. (2010) and Henze et al. (2011). A nova 
explosion is thought to be caused by a thermonuclear runaway 
in the hydrogen layer, which is accumulated on the surface 
of the white dwarf via accretion from the outer envelope of a 
companion star. Mainly depending on the mass of the white 
dwarf and the accretion rate of the system there can also be a 
state where the accumulated matter on the white dwarf surface 
undergoes steady hydrogen burning without a nova explosion 
(see |Nomoto et al.|2007[ ). We detect six SSSs in the catalogue, 
which are not associated with known optical novae (No. 64, 
184,229,294,295,314). 

The catalogue contains 45 X-ray sources associated with optical 
novae, 12 of which have not previously been published and will 
be discussed in detail by Henze et al. 2013 (in preparation). 
One type of novae shows very short outbursts in the light curves 
(12 classified as OB sources), while others show a very steady 
X-ray outburst over several years. Due to these long outbursts 
sources No. 79, 100 and 123 are classified as HV sources, as 
they show the longest outbursts we observe for novae (e.g. Fig. 
10 111. No. 123 even shows a re-brightening in the 2011/2012 



monitoring observations (see Fig.pLO|i 
Fig. 



11 gives an overview of SSSs and novae classified as 



significantly time variable in the catalogue. It can be seen, that 
most novae are classified as OB sources. There are three OB 
SSS without optical nova association (No. 184, 295, 64), which 
might be X-rays from a nova missed in optical obsei-vations. 
No. 294 and 229 are HV SSSs and show a particularly low 
variability factor, which as described above could be caused by 
steady hydrogen burning on the surface of a white dwarf. 



Discussion of individual SSS: 

No. 64 is a SSS showing an SOB, which is not associated 
with a reported optical nova (see |Henze et al.| |2010). However, 
the nova outburst might have been missed in optical observa- 
tions. 

No. 81 and No. 150 are classical novae in a M31 globular 
clusters (see| Henze et al.|2009b|[20T3l l. 

No. 184 is not detected in a 29.07 ks ACIS-I observation (ObsID 
7064 on 2006-12-04), where it should be bright according to 
the HRC-I light curve (see No. 184 in Fig. E\, which shows 
a plateau during this time. This indicates a super-soft X-ray 
spectrum, since the Chandra ACIS detector has a lower sensi- 



tivity than HRC-I at lower energies (see Pietsch et al. 2007 1. 
For a 50 eV black-body spectrum, the conversion factor from 
HRC-I to ACIS-I count rates would be 0.081 according to 
WebPIMMS. Thus converting the HRC-I (ObsID 7285) count 
rate of 1.95 ■ 10" "* counts s"', we would still expect 4.5 photon 
counts during the ACIS-I observation. This means, that either 
the spectrum is even softer, than we estimate above, or the 
source was unexpectedly fainter during the time of the ACIS-I 
observation and then re-brightened. The HRC-I light curve is 
typical for a SSS phase of an optical nova several months after 
outburst, but no known optical nova correlates with its position. 
No. 205 is detected a t the position o f nova M31N 2002-08b 
(optical light curve see|Lee et al.|2012l. In two Ha images from 



the [Massey et al. (2006 
Ho- emission visible at t 



I survey on 2002-09-12, there is bright 
le source position, which is most likely 

associated with the nova. This is the first detection of X-rays 

from M31N 2002-08b (see Table[TO|. 

No. 229 is a super-soft source candidate, showing a 217s pulse 

period (see Trudolyubov & Priedhorsky 2008|. 



No. 294 is a HV and classified as SSS by Pietsch etal.] ( |2005b| l, 
but does not show a nova outburst and is detected throughout 
the thirteen years of monitoring. It also has a optical counterpart 
in an Ha (see Table Hll image of the LOGS. 
No. 295 is a super-soft source, which could be the counterpart 
of a nova explosion missed in the optical. It has a 865s period, 
which is discussed in Osborne et al. ( 2001|l. Its spectrum is 
similar to galactic super-soft sources (see Trudolyubov et al. 

2001) . 

No. 314 shows a faint SOB (see [Haberl et aL][2007| ). The 
source re-brightened in a Swift follow-up observation and has 
a super-soft spectrum according to Kong & Di Stefano ( 2007[ l. 
This may be emission from a missed optical nova. 



4.2. Supernova remnants (SNR) 

After a supernova explosion, a shock front propagates through 
the interstellar medium (ISM), heating it to temperatures where 
thermal X-rays are emitted. There can also be strong radio 
emission from synchrotron radiation and optical emission lines 
are seen, where the ISM is ionized. Depending on the age of 
the SNR, the extent can be large enough to be resolved as an 
extended source of X-rays, radio or optical light even at the 
distance of M31. The catalogue contains six sources, thought 
to be associated with supernova remnants from the |Sasaki| 
et al. ( 2012[ l catalogue (some of their sources, overlapping 
with the FOV are not contained in the catalogue due to lower 
sensitivity of the HRC-I detector (compared to XMM-Newton 
EPIC) toward the edges of the FOV). Aside from those, we 
classified one new SNR candidate (No. 272), which shows 
X-ray extent, radio emission and hardness ratios (obtained from 
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Fig. 12. Contour plots of the logarithmically scaled merged HRC-I images with dashed circular regions at the radio positions with 
approximate radio radii ( Fesen et al.|1999 Sjouwerman & Dickel 2001) of the known supernova remnants within 90" of the center 
of M31. The images are smoothed with a Gaussian function of 7 pixel FWHM and contour levels are at 0.53, 0.56, 0.6 and 0.7 
counts per pixel. North is up and East is left. The zoom is slightly different (for region size and SNR names see Tablel9]). The X-ray 
catalogue source regions are circles with ~ 90% encircled count fraction radius. 
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image we find 446 source counts in a visually determined region 
(see Table|9|, which coiTesponds to a significance of about 3.2cr. 
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Fig. 11. Variability plot showing the spread of luminosity (left y- 
axis) and count rate (right y-axis) for all OB or HV SSS or nova 
sources. The x-axis is labelled with the catalogue source number 
(No.). If the minimum count rate is derived from an upper limit 
of the light curve, the lower count rate is marked by a downward 
arrow. The blue diamonds mark the sources which are classified 
as OB sources in the catalogue. The black dots mark the sources, 
which are associated to optical novae. 



Stiele et al.|2011 1 consistent with thermal X-ray emission. No. 
304 is a confirmed SNR and is discussed as being associated 
with a low-mass X-ray binary by Wil liams et al.| ( |2005a| l. 
By visually inspecting the merged high resolution image and 
smoothing with a 7 pixel FWHM Gaussian function, we found 
extended X-ray sources at the positions of the three radio SNRs 
discussed in [Sjouwerman & Dickel| ( |200 1 ) (see Fig. 12 1. Of 
these remnants only SNR 101 has a counterpart in the source 
catalogue (No. 215). SNR 101 and SNR 95 have alrea d y been 
detected in X-rays and are discussed in Kong et al. (2003[l. 



However our deeper image of SNR 95 is much more convinc- 
ing. For SNR 80 ([B90]80 in Braun 1990) we report the first 
observation in X-rays (Fig. 12 leftmost image), close to another 
bright X-ray source. In the 1.02 Ms merged high resolution 



In the right-most image of Fig. 12 we see a very dim X- 
ray emission at the position of SN 1885A ( [Fesen et al.||1999] l, 
which is not included in the catalogue. Source properties were 
extracted from a circular region with a radius of 0.7" (see 
estimated extent in [Sjouwerman & Dickel] |2001[ ) from the 
merged high resolution image at the position of SN 1885 A. 
We obtain 28 counts and a significance of 2.6cr, which is not 
sufficient to state a clear detection. To the south-east of SN 
1885 A, lies a transient X-ray source, which we found not to be 
associated with the supernova in accordance to Kaaret ( 2002[ l, 
since the positional eiTors plus the maximum extent of the 
source are clearly separated. Table l9] shows the circular regions 
we estimated for each SNR and the counts and significance we 
obtained for each of these regions using the CIAO tool aprates 
to compute source statistics. 



4.3. Foreground stars (fgStar) and background galaxies 
(AGN, Gal, GCI) 

There are sources in the catalogue, which are not located in 
M31. There are X-ray-bright foreground stars in our own 
galaxy, which we identified by correlating the source positions 
with known foreground stars (fgStar) in other X-ray catalogues 
(see section |3.1| l and visually inspecting images from the 
LGGS for clearly visible foreground stars. The advantage of 
the catalogue for these identifications is, that they get more 
precise for most sources, due to the good positional accuracy 
of Chandra HRC-I, and the long exposure available to us. The 
catalogue contains four confirmed fgStars and seven fgStar 
candidates, which are statistically distributed across the FOV as 
expected. For No. 85, 262, 289 and 291 we give the colours of 
the optical counterparts in Table |4] 

The luminosities of fgStars are eiToneously overestimated in the 
catalogue, since they are much closer to us than M 3 1 and we 
used the distance to M 31 for the calculation of all luminosities. 
Active galactic nuclei (AGN) in the background should show 
no X-ray extent in the images, have a hard spectrum and 
show moderate variability (see [Barnard et al.|2012b[ ). They are 
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Table 9. Source statistics for SNRs within 90" of the center of M 3 1, extracted from the 1 .02 Ms merged high resolution image. 



RA(J2000) [h:m:s] Dec(J2000) [d:m:s] Counts* Significance* Diameter ["] ' 



SNR ID " 



SNR80 


00:42:40.37 


+41:15:52.2 


446 


3.2 


7.8 


SNR 95 


00:42:47.81 


+41:15:26.3 


731 


4.5 


8.7 


SNR 101 


00:42:50.44 


+41:15:56.4 


856 


6.7 


7.9 


SNR 1885 A'' 


00:42:42.89 


+41:16:05.0 


28 


2.6 


1.4 



Notes. *°' Taken from Sjouwerman & Dickel J200lf . '*' Calculated with CIAO tool aprates. *'^' Visually estimated from smoothed merge image. 
*''' Position taken from Fesen et al.li 1999 1. 
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Fig. 13. Light curve of No. 99 during HRC-I observation 13178, 
showing an X-ray flare at the end of the observation. The plot 
displays counts over time in seconds from the start of the obser- 
vation, with a binning of 100s. 



super-massive black holes at the centre of a galaxy, accreting 
matter. The real X-ray luminosity is underestimated in the 
catalogue, since those sources are much further away from us 
than M31. AGN can also show strong radio emission. The 
radio emission in AGN is expected to arise from synchrotron 
radiation connected to acceleration processes in the jet of the 
super massive black hole (see Cavagnolo et al. 2010|). The 
catalogue only contains one confirmed AGN and one candidate, 
which we adopted from previous X-ray catalogues (section 3. 1 1. 



They show very faint light curves and thus we cannot draw any 
conclusions about variability. 

We did not find any identifications of X-ray sources with normal 
galaxies (Gal) or galaxy clusters (GCl) in the background nor 
with CVs or other known sources in the foreground of the M 3 1 
central field. 



Discussion 
objects: 



of individual foreground and background 



No. 44 and 45 show radio emission (see B90: |Braun|1990| l, and 

a harder spectrum ( jStiele et al.|2011 1, which would suggest they 

are AGN candidates. 

No. 85 correlates with the position of a confirmed foreground 

star from the RBCv5 (see'G alleti et al.|2007 1 and the light curve 

shows one significant outburst, which is caused by a short flare 

in the fight curve of HRC-I ObsID 10882. 

No. 99 is a fgStar candidate with time variability due to flaring 

activity in HRC-I ObsID 13178 (see observation light curve in 

Fig.[n]). 

No. 207 shows an optical counterpart (Table HI) and radio 
emission (source [B90]96). From the two extended optical 
sources, the radio emission and the soft X-ray spectrum at the 



position, the source could be a galaxy merger in the background, 
where the X-ray emission is caused by gas of the ISM, heated 
during the collision. 



4.4. M31 core 

For the largely overlapping central sources of M 31 we propose 
a new distribution of discrete X-ray sources in the catalogue (see 
Fig. [3}. 

No. 164 and 165 are at the position of the double core of M 31 
(see Ga rcia et al.||2010 l. They show variability very similar to 
the XRB light curves and are classified as HV sources in the 
catalogue. In Hubble Space Telescope images of the M31 core 
(M31*), a double structure, corresponding to the two X-ray 
source positions of No. 164 and 165 is visible as the most lu- 
minous region around the dynamical centre of the galaxy. This 
circum-nuclear environment is discussed as a region of very high 
star density around a central black hole by [ Bender et al.| ( [2005[ l 
and |Li et al. (2009). Li et al. ( 201 1| ) recently discussed No. 164 as 
the X-ray source of the central super-massive black hole (M 31*) 
and give a detailed Chandra light curve with an outburst in X- 
rays. The maximum of this outburst has only been observed with 
Chandra ACIS and is thus not visible in the HRC-I light curve 
of No. 164. 

No. 167 is a very soft X-ray source to the south of the core of 
M31 (see |Garciaetal.|2000| ). 

No. 170 and 172 are close to 167 and show cross talk from this 
source, so that we cannot classify them separately. 



5. Summary and conclusions 

We present the deepest X-ray catalogue of the M 3 1 central field 
available to date, which contains 318 discrete X-ray sources with 
detailed variability information (spanning 13 years). From the 
HRC-I data analysis, we are able to provide a crude classifica- 
tion of X-ray binary (XRB) candidates based on their long-term 
time variability. The catalogue contains light curves for 97 XRB 
candidates in the field of M 31 and 18 XRB candidates in globu- 
lar clusters (GIC). 

This classification of XRBs from time variability is the major 
step forward of the HRC-I analysis compared to previous stud- 
ies. Many of the XRB candidates, had only been classified as 
< hard > X-ray sources in one of the most recent X-ray cata- 
logues of M 31 byjStieIeetar| ([2011]). 

36 % of all sources in the catalogue are classified as XRB candi- 
dates with hard X-ray spectiTim (X-rays from optical novae make 
up another 14 % of the catalogue sources). Nine XRB candi- 
date sources show recurrent outbursts. The analysis may suggest, 
that outburst sources are less frequent in GlCs than in the field 
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of M31. We detected one SNR in the central region of M31, 
which is seen in X-rays for the first time and discovered one pre- 
viously unknown SNR candidate. We present some promising 
optical counterparts to the X-ray sources, which could be spec- 
trally identified with future optical follow-up observations. 
We observed 61 new sources, which have not yet been pub- 
lished in a dedicated M31 catalogue. 17 of those sources have 
never been detected before, 12 sources are new X-ray sources 
associated with optical novae, 19 sources have only been de- 
tected in the Chandra source catalogue (see Evans et ar]|2010| l 
and 13 sources have only been discussed in separate papers or 
Astronomer's Telegrams, but are not included in previous cata- 
logues. 

The catalogue contains a large sample of XRBs, with variabil- 
ity behaviour and can thus provide important input for a better 
understanding of XRBs and the X-ray source population of the 
M 3 1 central field in general. 
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Table 10. New X-ray 


sources not previously published in dedicated M 31 catalogues (section 3.2 . 
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Dec [d:m:s] 
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CSC name Ref. '' 


Comment 


11 


00:41:57.74 


+41 
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48.56 


3.56E+38 


26.09 
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+41 


10:30.8 


10.26 


1.23E+37 


4.96 


J004208.3+411029 
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00:42:11.25 


+41 
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13.86 


6.20E+37 
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2 
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00:42:17.29 
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315.73 


2.92E+38 


42.67 


J004217.3+411537 
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47 


00:42:21.29 
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48.01 


3.94E+37 
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AP 


65 
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15:47.6 
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8.84E+36 


7.29 


J004229. 1+41 1547 
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67 


00:42:29.10 


+41 


13:48.7 
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AP 
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J004231.7+411454 






82 
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30.70 
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9.28 
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00:42:33.26 
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5.21 


1.27E+35 
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88 
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1.93 






M31N2009-05b'' 


94 
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15:54.8 


7.00 


9.66E+35 


1.96 


J004236.5+411554 






103 
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00:42:38.67 


+41 
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1.90 








220 
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+41 
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00:42:53.17 
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2.99 






M31N2011-10d'' 


244 


00:42:55.99 
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3.59 
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Fig. 14. HRC-I merged high resolution image (0.13" pixels) with source positions and source No. from the catalogue. The region 
size is the 90% encircled count fraction PSF. 
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Fig. 15. HRC-I merged middle resolution image (0.52" pixels) with source positions and source No. from the catalogue. The region 
size is the 90% encircled count fraction PSF. 
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Fig. 16. HRC-I merged low resolution image (2.08" pixels) with source positions and source No. from the catalogue. The image 
shows the overlap of the whole detector at different roll angles. The region size is the 90% encircled count fraction PSF. 
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Table A.l. Chandra HRC-I observations used in this study. 



ObsID 


Exposure[s] " 


RA(J2000) [deg] 


Dec(J2000) [deg] 


Roll angle [deg] 


MJD* 


UT^ 


00267 


1268 


10.68195 


41.26369 


264.83 


51512.77 


1999-11-30.77 


00268 


5177 


10.68383 


41.26435 


282.14 


51535.10 


1999-12-23.10 


00269 


1204 


10.68627 


41.26340 


298.54 


51562.91 


2000-01-19.91 


00271 


2463 


10.68926 


41.26492 


331.70 


51611.76 


2000-03-08.76 


00272 


1202 


10.68767 


41.27255 


82.55 


51690.70 


2000-05-26.70 


00273 


1190 


10.68734 


41.27320 


101.67 


51716.20 


2000-06-21.20 


00275 


1188 


10.68291 


41.27314 


134.00 


51774.43 


2000-08-18.43 


00276 


1182 


10.68046 


41.27222 


152.42 


51798.86 


2000-09-11.86 


00277 


1187 


10.67952 


41.26924 


192.16 


51829.18 


2000-10-12.18 


00278 


1182 


10.68094 


41.26537 


250.96 


51865.19 


2000-11-17.19 


01570 


1183 


10.68592 


41.27310 


94.69 


52070.93 


2001-06-10.93 


01912 


46731 


10.67112 


41.26732 


226.83 


52213.99 


2001-10-31.99 


02904 


1181 


10.68094 


41.26564 


253.90 


52232.88 


2001-11-19.88 


02905 


1094 


10.68620 


41.26497 


296.74 


52290.20 


2002-01-16.20 


02906 


1187 


10.68861 


41.27227 


88.37 


52427.90 


2002-06-02.90 


05925 


46288 


10.68295 


41.26512 


270.75 


53345.76 


2004-12-06.76 


06177 


20037 


10.68427 


41.26494 


285.27 


53366.31 


2004-12-27.31 


05926 


28265 


10.68433 


41.26482 


285.63 


53366.84 


2004-12-27.84 


06202 


18045 


10.68601 


41.26498 


303.50 


53398.06 


2005-01-28.06 


05927 


26999 


10.68605 


41.26490 


303.50 


53398.83 


2005-01-28.83 


05928 


44855 


10.68742 


41.26528 


319.07 


53422.68 


2005-02-21.68 


07283 


19941 


10.68668 


41.27304 


90.26 


53891.29 


2006-06-05.29 


07284 


20002 


10.67974 


41.27133 


173.77 


54008.89 


2006-09-30.89 


07285 


18516 


10.68043 


41.26645 


245.26 


54052.30 


2006-11-13.30 


07286 


18707 


10.68829 


41.26613 


333.63 


54170.61 


2007-03-11.61 


08526 


19944 


10.67992 


41.26685 


236.71 


54411.64 


2007-11-07.64 


08527 


19980 


10.68098 


41.26611 


250.73 


54421.76 


2007-11-17.76 


08528 


19975 


10.68183 


41.26549 


262.91 


54432.79 


2007-11-28.79 


08529 


18922 


10.68262 


41.26536 


270.77 


54441.57 


2007-12-07.57 


08530 


19866 


10.68330 


41.26507 


278.33 


54451.49 


2007-12-17.49 


09825 


20212 


10.67993 


41.26690 


234.91 


54778.34 


2008-11-08.34 


09826 


19918 


10.68095 


41.26603 


250.89 


54787.14 


2008-11-17.14 


09827 


19962 


10.68197 


41.26556 


263.10 


54798.24 


2008-11-28.24 


09828 


19967 


10.68257 


41.26534 


271.22 


54807.41 


2008-12-07.41 


09829 


10075 


10.68351 


41.26532 


279.41 


54818.02 


2008-12-18.02 


10838 


10041 


10.68331 


41.26521 


279.41 


54818.49 


2008-12-18.49 


10683 


19897 


10.68687 


41.26546 


315.66 


54878.90 


2009-02-16.90 


10684 


18697 


10.68756 


41.26565 


322.39 


54888.17 


2009-02-26.17 


10882 


18846 


10.68010 


41.26661 


239.91 


55142.23 


2009-11-07.23 


10883 


18326 


10.68083 


41.26613 


249.41 


55151.24 


2009-11-16.24 


10884 


18358 


10.68186 


41.26558 


262.25 


55162.63 


2009-11-27.63 


10885 


18267 


10.68271 


41.26531 


272.30 


55173.94 


2009-12-08.94 


10886 


18343 


10.68341 


41.26513 


278.92 


55182.90 


2009-12-17.90 


11808 


17114 


10.68681 


41.26536 


314.83 


55242.86 


2010-02-15.86 


11809 


18415 


10.68747 


41.26558 


322.27 


55253.27 


2010-02-26.27 


12110 


19954 


10.68059 


41.26639 


246.42 


55514.17 


2010-11-14.17 


12111 


19866 


10.68144 


41.26587 


257.32 


55523.18 


2010-11-23.18 


12112 


19944 


10.68230 


41.26531 


267.64 


55533.66 


2010-12-03.66 


12113 


18974 


10.68299 


41.26516 


274.87 


55542.56 


2010-12-12.56 


12114 


19977 


10.68367 


41.26520 


281.67 


55552.18 


2010-12-22.18 


13178 


17453 


10.68688 


41.26540 


315.50 


55609.15 


2011-02-17.15 


13179 


17482 


10.68753 


41.26564 


322.79 


55619.25 


2011-02-27.25 


13180 


17023 


10.68819 


41.26608 


332.15 


55630.12 


2011-03-10.12 


13227 


19987 


10.68029 


41.26644 


243.25 


55877.10 


2011-11-12.10 


13228 


19003 


10.68122 


41.26587 


254.80 


55886.22 


2011-11-21.22 


13229 


19567 


10.68213 


41.26552 


264.98 


55895.98 


2011-11-30.98 


13230 


18932 


10.68288 


41.26521 


273.93 


55906.56 


2011-12-11.56 


13231 


19487 


10.68355 


41.26516 


280.22 


55915.33 


2011-12-20.33 


13278 


18973 


10.68683 


41.26548 


315.75 


55974.76 


2012-02-17.76 


13279 


18825 


10.68765 


41.26564 


323.40 


55985.26 


2012-02-28.26 


13280 


19285 


10.68859 


41.26629 


335.95 


55999.21 


2012-03-13.21 


13281 


18907 


10.68685 


41.27307 


87.90 


56079.90 


2012-06-01.90 



Notes. *''•' Exposure after good time interval correction. **' Modified Julian Date at the start of the observation. *^' Universal Time at the start of the 
observation. 
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Fig.B.l. Light curve plots showing X-ray luminosity (left axis, 0.2-lOkeV) and photon count rate (right axis) over modified Julian 
date for all catalogue sources. Fluxes with \cr errors are given as black dots and 3cr upper limits as red triangles. 
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Fig.B.l. continued. 
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Fig.B.l. continued. 
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Fig.B.l. continued. 
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Fig.B.l. continued. 
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Fig.B.l. continued. 



